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Abstract 

By the end of the current decade, there are plans to deploy several geostationary 
Earth orbit (GEO] satellite missions for atmospheric composition over North 
America, East Asia and Europe with additional missions proposed. Together, these 
present the possibility of a constellation of geostationary platforms to achieve 
continuous time-resolved high-density observations over continental domains for 
mapping pollutant sources and variability at diurnal and local scales. In this paper, 
we use a novel approach to sample a very high global resolution model (GEOS-5 at 7 
km horizontal resolution] to produce a dataset of synthetic carbon monoxide 
pollution observations representative of those potentially obtainable from a GEO 
satellite constellation with predicted measurement sensitivities based on current 
remote sensing capabilities. Part 1 of this study focuses on the production of 
simulated synthetic measurements for air quality OSSEs (Observing System 
Simulation Experiments]. We simulate carbon monoxide nadir retrievals using a 
technique that provides realistic measurements with very low computational cost. 
We discuss the sampling methodology: the projection of footprints and areas of 
regard for geostationary geometries over each of the North America, East Asia and 
Europe regions; the regression method to simulate measurement sensitivity; and 
the measurement error simulation. A detailed analysis of the simulated observation 
sensitivity is performed, and limitations of the method are discussed. We also 
describe impacts from clouds, showing that the efficiency of an instrument making 
atmospheric composition measurements on a geostationary platform is dependent 
on the dominant weather regime over a given region and the pixel size resolution. 
These results demonstrate the viability of the "instrument simulator” step for an 
OSSE to assess the performance of a constellation of geostationary satellites for air 
quality measurements. We describe the OSSE results in a follow up paper (Part 2 of 
this study]. 
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1. Introduction 

Current satellite observations of tropospheric composition made from low 
Earth orbit (LEO) satellites provide at best one or two measurements each day at 
any given location. Coverage is quasi-global but sparse, often with large 
uncertainties in individual measurements that limit examination of local and 
regional atmospheric composition over short time periods. This has hindered the 
operational uptake of these data for monitoring air quality and population exposure, 
and for initializing and evaluating chemical weather forecasts. 

By the end of the current decade, there are planned geostationary Earth orbit 
(GEO) satellite missions for atmospheric composition over North America, East Asia 
and Europe, with additional missions proposed (CEOS, 2011). Together, these 
present the possibility of a constellation of GEO platforms to achieve continuous 
time-resolved high-density observations over continental domains for mapping 
pollutant sources and variability. The GEO geometry provides a continuous view of 
the part of the Earth that is below the satellite, enabling measurements many times 
per day that help capture the diurnal evolution of emission sources, tropospheric 
chemistry and pollution transport. 

The planned GEO missions include the EVI-l/TEMPO (Tropospheric 
Emission: Monitoring of Pollution, Zoogman et al., 2014b) over USA, Sentinel 4/IRS 
over Europe and GEMS over Asia. Each mission has a different primary objective, 
but they share the common goal of monitoring pollutants for atmospheric 
composition and air quality and will have a common measurement capability for 
ozone (O 3 ), nitrogen dioxide (NO 2 ), sulfur dioxide (SO 2 ), formaldhyde (HCHO) and 
aerosols, utilizing radiances in the ultraviolet-visible (UV-Vis) spectrum. Planned 
GEO observations of infrared active trace gases of relevance to air quality are 
currently limited to total column carbon monoxide (CO) observations from the 
European IRS instrument, which is originally not driven by atmospheric 
composition applications. However, other IR measurements that could play a part in 
the GEO constellation are being proposed as part of the NASA Decadal Survey GEO- 
CAPE (GEOstationary Coastal and Air Pollution Events) mission, such as the EVI-3 
CHRONOS mission ( https: //www2.acd.ucar.edu/chronos ) that would measure CO 
and methane (CH 4 ) using heritage from the Terra/MOPITT (Measurement of 
Pollution in The Troposphere) instrument. Given the effect of nearby emissions and 
transported pollution on local air quality, MOPITT-like CO observations are a good 
candidate for air quality measurements on a GEO platform because the unique 
sensitivity of this platform to pollution in the boundary layer, as well as in the free 
troposphere, allows both vertical and horizontal tracking of pollution transport. 

Carbon monoxide is a primary pollutant and plays an important role in 
tropospheric chemistry and its sources are both natural and anthropogenic. There 
are two main processes of CO production: incomplete combustion (e.g., industrial 
and urban fossil/bio fuel burning, wildfires and biomass burning); and natural 
chemical production from hydrocarbon oxidation. As an O3 precursor, CO is also 
important in determining the tropospheric O3 budget. The principal CO sink is the 
oxidation by hydroxyl (OH) radicals, giving an average CO lifetime of about two 
months dependent on season. With these characteristics, CO serves as a tracer of 
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pollution emissions and transport, and as a proxy for emissions and distributions of 
other species co-emitted with CO but not easily measured. Taken together, 
observations of the full suite of UV-Vis and IR trace gases and aerosols could provide 
the high spatio-temporal resolution continental-scale observations of lower- 
tropospheric pollution needed to monitor, forecast, and manage air quality on a 
daily basis (Edwards et al., 2009; Lahoz et al., 2012; Bowman et al. 2013). 

Previous GEO observation simulation studies for air quality have assessed 
the potential capabilities of instruments covering the above three continental 
regions separately. Edwards et al. (2009) and Zoogman et al. (2011, 2014ab) 
consider the CONUS (continental US) region and demonstrate the feasibility of using 
observing system simulation experiment (OSSE) studies to help define quantitative 
trace gas measurement requirements in different spectral regions for satellite 
missions and to evaluate the expected performance of proposed observing 
strategies to test the ability of GEO satellite measurements of ozone (O 3 ) and CO. 
Claeyman et al. (2011) and Sellitto et al. (2013b) cover the European region and 
describe the capabilities of a concept nadir thermal infrared sensor proposed for 
deployment onboard a GEO platform to monitor O 3 and CO for air quality purposes 
(MAGEAQ: Monitoring the Atmosphere from Geostationary orbit for European Air 
Quality). Lastly, Zoogman et al. (2014a) assimilate concurrent ozone and CO 
observations and show that geostationary measurement of CO provides significant 
benefit for monitoring ozone. 

The goal of this study is to evaluate the impact of a future GEO constellation 
on global chemical weather by using the observing system simulation experiment 
(OSSE) technique. Here we primarily consider CO as a good chemical tracer for 
evaluating the impact of a GEO constellation of observations. As described by 
Edwards et al. (2009), chemical OSSEs provide a way of expanding case-specific 
sensitivity studies to assess the impact of future measurements systems. A chemical 
OSSE is composed of several elements (fig. 1 - see also, Timmermans et al., 2014). A 
nature run (NR) (1) represents the atmospheric true state. A complete OSSE needs 
an observation simulator (2) to sample the nature run to produce synthetic 
observations (3). The synthetic observations are then assimilated using a data 
assimilation system (4) into a second atmospheric model, the control run (CR) (5). 
This produces the assimilation run (AR) (6). The impact of concept instrument 
measurements on constraining the modeled state of the atmosphere can then be 
evaluated and assessed (7) by comparing the NR, CR and AR (1, 5 and 6). We 
describe this study in two parts. In the present paper (Part 1) we focus on the NR 
(1), observation simulator (2) and synthetic observations (3). A follow-up article 
(Part 2) will focus on assimilating the simulated measurements and assessing the 
synergies between the different instruments of the constellation by simulating data- 
denial case studies (elements 4 to 7 in fig. 1). This study presents for the first time a 
global GEO constellation OSSE for CO. 

According to Rodgers (2000), within the remote sensing optimal estimation 
framework one can represent the sensitivity of the retrieved trace gas profile from a 
satellite measurement to the true state of the atmosphere by the averaging kernel 
(AK) function. For accurate observation simulations in an OSSE, we need a full 
radiative transfer model for radiances and their Jacobians (which represent the 
sensitivity of the radiance to the true atmospheric state) to compute the AKs for 
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each atmospheric and surface scene. Since this presents a significant computational 
burden, practical implementations of OSSEs for air quality to date have used 
approximated observation simulators. Some have used specified constant AKs 
(Edwards et al., 2009; Zoogman et al., 2011), or have simplified the AK variability by 
considering only a few scene types (Claeyman et al., 2011). Sellitto et al. (2013a) 
showed that the use of no, or limited scene dependent AK, parameterizations could 
significantly misrepresent the sensitivity of an observing system. Sellitto et al. 
(2013a) also recommend using comprehensive scene-dependent approximations of 
the AKs in cases where the computational cost of a full radiative transfer model is 
too expensive to perform an OSSE study (for example, for a GEO constellation). 
Worden et al. (2013) address this issue by using a multiple regression analysis of 
real satellite observations to estimate scene-dependent averaging kernels, thus 
avoiding the use of a full radiative transfer model. This method allows the fast 
computation of scene-dependent AKs, and the processing of a very large dataset of 
synthetic observations in a short amount of time. 

Due to the constraint from the NR space and time resolution, approximations 
made to the instrument sampling and horizontal resolution cannot provide 
information at a higher resolution than the nature run (Edwards et al., 2009). One 
should use high space and time resolution NRs to simulate high instrument space 
and time sampling. The planned missions mentioned above would provide less than 
10 km spatial resolution at about every 1 hour. Sellitto et al. (2013b) also 
approximated the observation simulation by not discarding the cloud-contaminated 
measurements, thus leading to a possible overestimation of the GEO instrument 
potential to monitor tropospheric O 3 and pollution features in general. One should 
account for cloud contamination by testing scenarios with variable instrument 
sampling and resolution. 

In this paper we use the multiple regression analysis of Worden et al. (2013) 
to produce a very large data set representing a GEO constellation of synthetic 
observations for air quality. In section 2, we describe the very high resolution NR 
from the Goddard Earth Observing System Model version 5 (GEOS-5) at 7 km 
horizontal resolution. Section 3 describes the sampling methodology with details on 
the geostationary projection to the surface of the earth, and the multi linear 
regression method with its limitations for predicting averaging kernels and 
estimated observation errors. Section 4 investigates the impacts of clouds on the 
GEO constellation. The effect of horizontal resolution and sampling is discussed. 
Section 5 presents the measurements and a detailed analysis of the simulated 
observation sensitivity (e.g., averaging kernel variability). Section 6 gives a 
summary, conclusions and perspectives. 

2. The nature run 

The Goddard Earth Observing System Model, Version 5 (GEOS-5, Rienecker et 
al., 2008) is used to provide the NR. The GEOS-5 atmospheric model is a weather- 
and-climate model used for atmospheric analyses, weather forecasts, uncoupled and 
coupled climate simulations and predictions, and for coupled chemistry-climate 
simulations. The NR used for this study covers a 2-year global, non-hydrostatic 
mesoscale simulation for the period 2005-2006. In addition to standard 
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meteorological parameters [wind, temperature, moisture, surface pressure), this 
simulation includes 15 aerosols tracers [dust, sea salt, sulfate, black and organic 
carbon), and O 3 , CO and carbon dioxide [CO 2 ) trace gases. 

The model simulation is driven by prescribed sea-surface temperature and 
sea-ice derived at a horizontal resolution of 0.25 degrees. Biomass burning 
emissions of organic carbon, sulfate, CO and CO 2 are obtained from the Quick Fire 
Emissions Dataset [QFED) version 2.4-r6. The basis of the QFED is the fire radiative 
power [top-down) approach, and it draws on the cloud correction method used in 
the Global Fire Assimilation System [GFAS; Kaiser et al. 2012). Anthropogenic 
emissions of carbon species and aerosols are largely taken from the Emissions 
Database for Global Atmospheric Research [EDGAR; Olivier et al., 1994), which are 
provided annually at a resolution of 0.1 degrees. For CO and CO 2 , EDGAR v4.2 
emissions from 2005 through 2007 were used. For organic and black carbon 
aerosols species, Hemispheric Transport of Air Pollution (HTAP) emissions were 
used. 

Outputs at 30-minute intervals have been produced at a resolution of 0.0625 
degrees [~7 km) using a cubed-sphere horizontal grid with 72 vertical levels, 
extending from the surface up to 0.01 hPa [ — 85 km). All details and references 
concerning nature run file specifications, meteorology, chemistry and emissions can 
be found in the NR description documents at: 
http://gmao.gsfc.nasa.gov/proiects/G5NR/ 

In this study we focus on July 2006. Figure 2 shows the CO total column 
provided by the NR for 15 July 2006 at 15:00 UT. This map shows the ability of the 
NR to represent the high variability of CO fields at a global scale. We display typical 
and expected CO values: very high values [above 4.10 18 molecules/cm 2 ) over central 
Africa due to biomass burning; high values [around 3.10 18 molecules/cm 2 ) over 
dense populated areas due to anthropogenic emissions. The NR total columns of CO 
also clearly show long-range transport patterns of CO from anthropogenic and 
biomass-burning sources across the oceans of the Northern Hemisphere [NH) and 
Southern Hemisphere [SH), respectively. 

Figure 3 shows the July 2006 average of surface CO values over the three 
regions of interest [North America - CONUS, Europe and Eastern Asia). The NR 
shows realistic horizontal CO variability due to the very high space and time 
resolutions of the simulations. Emissions from cities from small to large size are 
clearly identifiable. Transport infrastructure such as roads [eastern US in figure 3.a) 
and ship routes [China sea in figure 3.c) are also visible. In this study we use the NR 
model output variables, both the chemical parameters [CO quantities) and the 
meteorological parameters [not shown), to predict averaging kernels for simulated 
observations in the GEO constellation. This is done for each of the CONUS, Europe 
and Eastern Asia regions of interest. 

3. Sampling methodology 

3.1. Geometry of measurements 

We constructed three GEO instrument simulators over the three regions of 
interest defined immediately above using the methodology described in Worden et 
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al. (2013). Footprints of the instruments are defined as a GEO projection on the 
globe. We defined x (along the parallel from the sub-satellite point) and y (along the 
meridian from the sub-satellite point) at regularly spaced scanning angles (in 
degrees). The GEO projection consists of projecting these angles from the GEO 
platform to the surface of the earth to obtain the corresponding longitudes and 
latitudes of the footprints. We have the following relationship between viewing 
angles at the satellite location and latitude, longitude position on the earth surface: 

Ion = tan _1 (s 1 /s 2 ) + subjon (1) 

lat = tan -1 (p 2 (s 3 /s xy )j (2) 

where subjon is the sub-satellite point longitude and: 

S 1 = Pi — s n cos x COS y 

s 2 = s n sin x cos y 
s 3 = —s n siny 



s d = >/(Pi cos x cosy ) 2 — ((cosy ) 2 + p 2 (sin y) 2 )p 3 

p 1 cos x cos y — s d 

s — 

n (cos y) 2 + p 2 (sin y) 2 

pi = 42164 km, the altitude of a GEO platform from the center of the earth 

P2 = 1.006803 is the ratio of the earth radius at the equator and at the pole (p2 = 

req/rpo). 

P 3 = pi 2 - r eq 2 

These equations follow from the methods provided in the technical 
document EUMETSAT (2011) and sketch of figure 4.d should be consulted to 
understand the above formulas. Projecting the regularly spaced instrument viewing 
angles onto the surface of the earth (figure 4.b) results in GEO instrument footprints 
with non-regular latitude and longitude spacing. GEO instruments then have a non- 
uniform horizontal resolution: the footprint density per surface area decreases as 
the measurements go outward from the sub-satellite point (figure 4.c). The GEO- 
CAPE concept mission (Fishman et al., 2012) requires hourly measurements with a 
spatial resolution in the order of 5 to 10 km and a measurement domain of at least 
5000 km. Table 1 gives an overview of the characteristics of the three instruments 
that we call hereafter GEO-US (over CONUS), GEO-EU (over Europe) and GEO-AS 
(over Eastern Asia). We set the scanning angles of the three instruments to have a 
horizontal resolution under 10 km (0.1°) in the approximate middle of the 
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measurement domain (i.e., sub-longitude and the mean of latitudes at the sub- 
longitude]. Figure 4.a. shows the measurement domains of the GEO constellation. 
Areas of coverage have different shapes due to the latitudinal extent of continents; 
GEO-EU has more of a latitudinal extent compared to GEO-US, which has to cover a 
wider longitude range. GEO-AS has been designed as a compromise solution 
between measurements over Chinese mega-cities and measurements over Korea 
and Japan. 

3.2 Carbon monoxide instrument simulator 

In this study, we assume characteristics of the CO measurements of the 
troposphere similar to those of the Terra/MOPITT (Measurement of Pollution in the 
Troposphere) instrument (Drummond, 1992). The last version of the retrieved CO 
product version 5 (Deeter et al., 2013) uses a multispectral approach utilizing near- 
visible infrared (NIR) solar backscatter signals at 2.3 microns and thermal infrared 
(TIR) emission signals from the Earth surface and atmosphere at 4.6 microns. This 
approach provides enhanced measurement sensitivity to near-surface CO 
concentrations and allows the possibility of retrieving CO profile information to 
separate CO in the planetary boundary layer and free troposphere (Worden et al., 
2010). This is a requirement for the GEO-CAPE concept mission (Fishman et al., 
2012) and it is generally desirable for air quality space remote observations to 
distinguish between local emissions and transported pollution at a given location 
(Lahoz et al., 2012). In the case of MOPITT, the combination of the TIR and NIR 
radiances significantly improves the sensitivity to the lower tropospheric CO for 
daytime land observations. For nighttime land and day/night ocean observations, 
only the TIR radiances contribute to the retrieval. 

The MOPITT-retrieved CO volume mixing ratios (VMRs) are on 10 pressure 
levels (surface, 900, 800, 700, 600, 500, 400, 300, 200, 100 hPa). Each retrieved 
level is representative of the layer content defined by the level value itself and the 
level above. The top most level extends from 100 hPa to 50 hPa. The retrieved CO 
profile y r can be related to the true atmospheric state y t with the following linear 
relationship: 

y r = y a + A{y t -y a ) + £ C3) 

In Eq. (3) y t is the true atmospheric CO profile state (in logio(VMR)) and y a is the a- 
priori state vector (in logio(VMR)) derived from a monthly mean climatological 
profile from the MOZART-4 (Model for Ozone and Related chemical Tracers, version 
4) chemical transport model (Emmons et al., 2010). The random errors (in 
logio(VMR)) is simulated using the retrieval noise, and A is the retrieval AK matrix 
(see section 3.3). (The y r retrieved profile obtained is then converted from 
logio(VMR) to VMR for the final data product). 

Figure 5 shows two representative MOPITT AKs. The sensitivity of the 
MOPITT instrument to near-surface CO varies according to different surface types 
and atmospheric conditions. The left panel of Fig. 5 shows a typical AK for a daytime 
measurement over land with enhanced sensitivity toward the surface. The right 
panel of Fig. 5 shows a typical AK for a TIR-only ocean or nighttime measurement 
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over land with low sensitivity in the lowermost troposphere. A useful quantity 
indicating the information content of a measurement is the degrees of freedom for 
signal [DFS], given by tr(A) [Rodgers, 2000]. Higher DFS values indicate more 
sensitivity of the retrieval to the true profile. 

To diagnose the sensitivity of the measurement to the lowest layers, DFS can 
be calculated over the three lowest levels [Surface to 700 hPa] as DFS 0 3 = £f =1 A a . 
In Figure 5 the DFS [DFSo, 3 ] is 1.9 [0.7] and 1.5 [0.2] for land-day and ocean-night 
measurements, respectively. We can see that MOPITT sensitivity toward the surface 
[DFSo, 3 ] is scene dependent. That is, it depends on various land and atmospheric 
parameters [i.e., nature of the surface and current state of the atmosphere at a given 
time] that control, among other things, the surface-atmospheric thermal contrast, 
i.e., the difference between the surface temperature and the atmospheric 
temperature profile. 

3.3 Simulated retrieval method. 

Worden et al. [2013] investigated the CO retrieval error resulting from the 
use of a single average AK in an observation simulator compared using the true 
retrieval AKs. They further developed a scene-dependent AK prediction tool capable 
of approximating the true AK with a significant reduction in retrieved CO error 
compared to using a single average AK. This AK prediction tool allows us to produce 
a large amount of simulated data over months in an efficient manner. One month of 
data for a GEO constellation [i.e,. around 200 million profiles] can be produced in 
less than 12 hours. 

The method of Worden et al. employs a multiple regression approach for 
deriving scene dependent AKs using predictors based on state parameters from the 
NR. The main parameters used are: CO concentration, temperature, specific 
humidity and pressure [see table 2]. The method is based on the computation of the 
singular value decomposition [SVD] of the AK matrix. Given an AK matrix A, we 
compute the SVD by means of: 

A = UAV r [4] 

where the columns of U and V are the left and right singular vectors respectively, 
and the elements of A (a diagonal matrix] are the singular values. Since the first two 
singular vectors account for 95% of the variability of MOPITT CO AKs on average 
and the first three singular vectors account for 99.995 %, the method retains the 
first three ranked singular vectors. For a complete description of the SVD technique, 
numerical examples and software used please refer to Worden et al., 2013. We then 
calculate the three first singular vectors and values using multiple regression. For 
example, 


Uij — Cij + 2fc=i Q-ijkXjk U) 

with dimensions of: / singular vectors, j pressure levels, and k predictors. The 
parameters are: c, a constant; a, regressions coefficients; and x, predictors. We used 
twelve predictors [Al=12] and have defined eleven different training sets [containing 
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the a coefficients] for the geographical regions of interest. Only a single training set 
can be used in the regression calculation. The predictors and training sets are listed 
in Table 2. Worden et al. (2013] selected the predictors based on their importance 
in the regression technique for parameterizing MOPITT forward model 
transmittances of Edwards et al. (1999]. The training sets are derived from a multi 
linear fit using real MOPITT observations. The training set period is the entire year 
2006. Once an AK matrix A is predicted, the simulated observation profile from the 
NR can be computed using the retrieval equation: 

y r = y a + My nr -y a ) + £ ( 6 ] 

with y NR , the NR profile sampled at the MOPITT vertical resolution, replacing the 
true state profile y t in equation 3. Because MOPITT retrieved values express a CO 
quantity over a pressure layer, we compute a weighted average using the pressure 
thickness of the GEOS-5 vertical CO levels mapped onto the MOPITT grid to produce 
Vnr- 

3.4. Training set method limitations 

In section 3.3 we applied the method described Worden et al. (2013), to 
reconstruct the averaging kernel matrix. In order to utilize the multi linear 
regression (equation 5), we need pre-calculated coefficients (ao,N) from a multi- 
regression fit derived from real MOPITT observations, that we call training sets 
given in Table 2. In some cases, mostly over the CONUS and Asian megacities, very 
high CO profile concentrations and total CO column amount values can extend 
beyond the boundary values of the data set used to build the training set and hence 
beyond the boundary values of the training set itself. Because of the near linear 
relationship between predictors and predicted AKs (equation 4 and 5), using 
predictors from the model with values that are outside the training set distributions 
may lead to unphysical averaging kernel values, e.g., strong negative values or 
values above unity. This is most likely the case for the CO predictors (CO profile and 
CO total column). In order to prevent predictors that are outside the training set 
range and not to discard a significant amount of simulated observation over 
polluted areas we reduce the CO profile predictor as follows. We calculate the mean 
(p) and standard deviation (a) of the CO profile training sets. If the predicted CO 
profile values (p) are above p+2 a, the new predictor (p’) is then calculated as 
follows: 

P' = (1-Y)(B + 2ct)+yp (7) 

where y is a weighting coefficient ranging between 0 and 1. Then the scaled CO 
profile is used to recalculate related CO predictors (CO column, Cos(0 S za)/logioCO(z) 
and dT(z) )/logioCO(z)). This allows the simulator to produce reasonable variability 
in measurement sensitivity while still including the high CO cases and without 
generating unphysical averaging kernels. Sensitivity tests during extreme pollution 
events have shown that using y>0.2 produces an unacceptably high frequency of 
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unrealistic averaging kernel functions. In order to have a robust observation 
simulator which does not produce unphysical averaging kernel values we use y=0.1. 

3.5. Simulated error method 

The regression method described above does not account for simulating 
measurement error (represented by the retrieval error covariance matrix] and 
retrieval noise. In order to simulate the error terms, we use the relationships 
between the AK matrix and the associated retrieval errors terms (Rodgers, 2000). 
The associated retrieval noise a is defined using the retrieval noise covariance 
matrix C n , derived from the retrieval error covariance matrix C x . Where a is the 
vector containing the square root of the diagonal elements of C n , and y err the 
vector containing the square root of the diagonal elements of C x . The retrieval error 
covariance matrix C x can be decomposed as the sum of two matrices (Deeter et al. 
2011 ): 

• A smoothing error covariance matrix C s that describes the expected error 
arising from differences between the true profile and retrieved profile, and 
which are due to the characteristics of the weighting functions and the 
influence of the a priori covariance matrix. 

• A retrieval noise error covariance matrix C n that quantifies the expected 
errors due to errors in the radiances. 

Then 

C x = C s + C n (8) 

with C s approximated using the a priori covariance matrix C a , as follows 

C s = (/ - A)C a (I - A) T (9) 
and C x directly calculable from C a and A 

C x = (J — A)C a (10) 

so that 

C n = C x (I-0-A) T ) (11) 

C n is mostly lower than C s but not negligible (see section 5.2 and figure 11). 
Relatively to C x , C n will increase if C s decrease (if A tends to be the identity /). We 
define C a as for the MOPITT v4 and v5 products (Deeter et al., 2010). The a priori 
covariance matrix C a incorporates the same variance value C 0 at all levels, with a 
constant correlation height P c over a pressure level p defining the off-diagonal 
elements. Thus, 
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Cat) = Coe-^-wJ ( 12 ) 


with P c =100 hPa and Co=(0.3 logioe) 2 . In order to simulate the random errors, we 
add a pseudo-random noise on each nature run sampled by a predicted AK: 


where ° denotes the Schur product and K (0, /) a matrix following a normal 
distribution of means equal 0 and standard deviation equal the identity matrix /. We 
also calculate the retrieval error profile as follows: 


Because the smoothing error C s mostly dominates on the error budget (equation 8), 
the impact of the random error e is low compared to the retrieval error profile and 
hence the accuracy of the retrievals are not significantly impacted. 

4. Impact of clouds 

Under cloudy conditions, the simplest approach for MOPITT-like 
measurements on a GEO platform would be simply to discard cloudy pixels and not 
perform retrievals. It is thus important to assess the impact of the cloud coverage on 
GEO measurements. In this study, a scene is considered clear when the interpolated 
cloud fraction from the NR is lower than 5% of a single footprint. This is the clear- 
sky condition used operationally with real MOPITT measurements. Cloud 
contaminated footprints with greater than 5% of cloud fraction would be discarded. 
Clouds properties are not used to predict AK variability. Figure 6 presents the ratio 
of cloud free pixels, over the month of July 2006 for the constellation. The ratio of 
cloud free pixels is the number of cloud free observations divided by the total 
number of possible observations (i.e., one per hour during one month] for a given 
pixel. Figure 9 gives an idea of instantaneous instrument coverage with a 5% cloud 
fraction threshold. The GEO-EU displays very few cloud-contaminated areas 
whereas the GEO-AS has very few cloud free areas. 

The cloud-free ratio geographical distribution shows differences between 
intra- and inter-continental regions. On average, GEO-EU has the highest ratio 
(60%) followed by GEO-US (40%) and GEO-AS (20%). Strong variations of the ratio 
are also observed for different weather regimes within each measurement domain. 
Mediterranean weather regimes such as western CONUS and the entire 
Mediterranean basin exhibit higher ratios, above 80%. Conversely, oceanic, 
subtropical and tropical regimes such as northern Europe, southern CONUS and 
southeastern Asia have lower ratios, below 20%. Over the GEO-AS field of view, 
Korea and Japan exhibit very low ratios around 10% due to East Asian monsoon 
effects that provide persistent convective cloud coverage. 

The value of the cloud free ratio depends on the spatial resolution of the 
observation (pixel size) and the cloud fraction threshold used. Figure 7 displays 
results of sensitivity tests on pixel size and cloud fraction threshold. We assume that 


£ = y r c ){ 2 ° J\f (o,/) 


(13) 



(14) 
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the lowest pixel size simulated is 7 km due to the model horizontal resolution. We 
can then increase the pixel size by averaging contiguous grid cells. It is shown here 
that with a given cloud fraction threshold, increasing the pixel size reduces the 
average cloud free ratio. We perform tests for varying cloud fraction thresholds to 
calibrate the assimilated data product. Variations of the cloud free ratio following 
variations in the cloud fraction threshold and the pixel size show the same patterns 
(but with a different range of values] for the three instruments of the constellation. 

To explain these patterns we display a specific case (figure 8] as an example 
of how the observed coverage changes with the two varying parameters. The case 
study presented shows two typical horizontal cloud structures: one of high 
granularity located over the eastern part of the plot, which is identified as 
convective structures, and the other of low granularity located on the northwest 
part of the plot, which is identified as a cold air front. Over low granularity areas, 
decreasing the cloud fraction threshold will not increase the cloud-contaminated 
area as much as it does over the high granularity areas. As an idealized example, one 
can imagine adding pixels around four single separated sparse pixels (a granular 
structure) and adding a pixel around a four-by-four pixel area (a non-granular 
structure). In the first case, there will be 8 pixels around each of the 4 original 
pixels, making a total of 32 additional pixels. In the second case 12 additional pixels 
will surround a 2 by 2 square. The increase in area will be larger with the granular 
structure than with the non-granular structure. 

In the more realistic case of our observation simulations, granularity can 
vary at different scales and at different times. We found that adjusting the cloud 
fraction threshold to 20% for a 42 km pixel size gives comparable statistics of cloud 
coverage as with the 5% threshold for a 7 km pixel size (see section 5.3 and figure 
12 ). 

5. Simulated GEO constellation measurements 
5.1. Simulated sensitivity analysis 

Figure 9 displays the maps of sampled Surface-700 hPa NR and retrieved 
partial columns and associated DFSo ,3 for the GEO constellation. Looking at DFS 0-3 
maps first shows that the observation simulator reproduces the variability of 
measurement sensitivity over the satellite measurement domains. The maps are 
snapshots during daytime, and show strong differences in DFS 0-3 between sea and 
land due to the different AK training sets used. The land training set simulates multi- 
spectral (TIR+NIR) retrieval AKs in contrast to the sea training set that simulates 
TIR-only retrievals. The DFSo ,3 values between land and sea surfaces are in 
agreement with figure 5: instrument sensitivities over land are generally higher 
than over sea, because of the availability of multi-spectral simulated retrievals. DFS 
variability over land, or over sea only, is also simulated using the multi-regression fit 
as described in section 3.3. To describe this variability, we will focus on the analysis 
over land. The most obvious variations of DFSo ,3 follow orography. The main reason 
is the reduction in the number of retrieved levels if surface pressure is lower than, 
e.g., 900 hPa. For a constant number of retrieved levels, the variation of the surface 
level layer thickness also plays a significant role (represented by the dP predictor; 
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see table 2), and a thinner surface layer will contribute less retrieval sensitivity. 
Variations of DFSo ,3 can also be correlated to the CO amount in the NR. This 
variability is represented with the CO total column and CO profile predictors. CO 
abundance is a strong predictor of sensitivity due to the use of logio(VMR) retrievals 
in MOPITT with corresponding weighting functions that have increasing magnitude 
for increasing VMR [Worden et al., 2013]. Finally, the temperature profile and 
thermal contrast (dT) play a significant role in the DFSo -3 variability, as expected for 
the TIR contribution in a multispectral instrument. While DFSo -3 depends more on 
predictors such as CO column and dP, all of the predictors in Table 2 add 
information to the regression fit, as tested in Worden et al. [2013). 

Figure 10 shows scatter plots of DFSo -3 versus the main DFS variability 
drivers, i.e., parameters mentioned above such as CO concentration, dP and dT. 
Night and day values are displayed [blue and red, respectively) showing the 
expected increase of sensitivity during day [simulating a multispectral retrieval) 
compared to night [simulating a TIR-only retrieval). For each region, using an 
alternation of day training sets and night training sets, designed to produce 
multispectral and TIR-only retrieval AKs, respectively, then simulates a diurnal cycle 
of sensitivity. Correlation of DFSo -3 with predictors gives an indication of which 
variables in the NR true state will drive measurement sensitivity. However, this is 
not a deterministic result since actual sensitivity depends on all the predictors, 
together with the distributions of those 

variables as compared to the training set distributions, indicated by the lines in 
Figure 10. Variations in the dependence on predictors can be seen by the different 
distributions in Figure 10 for CONUS, Europe and Asia. Over Asia and Europe, 
overall CO concentrations from NR show significantly lower as compared to the 
training set mean. For Asia, scatter plots do not show any clear dependence between 
DFSo -3 and CO concentrations. For Europe, the dependence is more marked during 
daytime. Lower CO predictor values compared to training set mean might lead to 
underestimation of DFSo -3 , however it fits a realistic range of values (from 0.25 to 
0.5). 

5.2 NR sampling and error budget 

The difference between the NR CO and the retrieved CO shows higher NR 
values than in the simulated retrievals (fig. 9). Retrieved values can be close to the 
NR if sensitivity (DFS) is high enough and/or the a priori CO profile is close enough 
to the NR. Cases with strong CO plumes in the NR can be identified in figure 9 over 
Asia (around 35°N and 115°E) and over Europe (around 5°E and 55°N). In the Asian 
case, the plume is very well detected in the synthetic retrieval, because over land 
GEO-AS has a DFSo ,3 above 0.5 and a priori profile concentrations close to the NR 
profile (not shown). In the European case, plumes are barely detected because over 
sea the GEO-EU has DFSo ,3 below 0.3 and the a priori profile concentrations are far 
from NR values. In general, retrieved CO concentrations are lower than the NR CO 
concentration because a priori values are lower than NR values. In certain cases (see 
fig. 9 for Asia around 110°E and 35°N), the opposite is observed; a priori 
concentrations are higher than the NR. The a priori profile, sampled from a lower 
resolution MOZART-4 climatology (see section 3.2) does not capture the specific NR 
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high-resolution features. Conversely, polluted areas are represented as relatively 
high CO over broad area, which can produce cases where y a is higher than y t . 

Figure 11 left panels show scatter plots of NR CO partial columns (Xt) versus 
retrieved CO partial columns (Xr) with night cases (blue) and day cases (red) over 
land. In general, night Xr values are farther from the Xt compared to the day Xr 
values. As explained in section 5.1 and in figure 10, DFSo ,3 values are lower during 
night than during day. Lower DFS will produce Xt values that are closer to the a 
priori (Xa). If Xa is far from Xt, the smoothing error (Xs) will increase with lower 
DFSo, 3 . Even if DFSo ,3 is high (around 0.7], Xs can be high if Xt is very far from Xa. In 
the case of GEO-US, values spread by 10-20 DU (Dobson Units] around the Xt=Xr 
axis, showing that Xa can be higher or lower than Xt. In the case of GEO-EU, the 
spread is lower because Xt is in general close to Xa. In the case of GEO-AS, Xr values 
are mostly lower than Xt values, showing that Xa is generally lower than Xt. 

Figure 11 right panel displays scatter plots of Xs (in % relative to Xr) versus 
the surface-700hPa partial column retrieval error (Xe). We see that Xe values are in 
the range expected from real MOPITT observations: between 15% and 30%. 
Following equation 8 and 14, diagonal values of C s should be lower or equal to 
diagonal values of C x and hence Xs should be lower or equal to Xe (if Xs is calculated 
as Xe). The condition is respected in most of the cases, but some Xs values are higher 
than Xe. Again, this happens when Xa is very distant from Xt, and due to the fact that 
Xa and Xt (i.e., y a and y t ) are not used in the calculation of the a priori covariance 
matrix (see section 3.5 and equation 9 and 12). The perfect estimate of C s would 
then be: 

C s = a~ A)(y a - y t )(y a ~ y t ) T (f ~ A) T (15) 

This can be estimated for this study since we are assuming the NR is the true state. 
However, for real observations it is not possible to estimate the actual smoothing 
covariance error matrix. Therefore, use of the method described in section 3.5 is 
more realistic, and will provide reasonable error estimates in most cases since Xs 
has generally lower values than Xe. 

5.3. Reduced resolution simulated observations 

In part II of this study, we will assimilate the simulated GEO-constellation 
into a global model. We will use the global chemistry - climate model CAM-Chem, 
including its full chemical scheme (Lamarque et al., 2012). State-of-the-art global 
atmospheric chemical models do not have high horizontal resolution. In this second 
part of the study, we use a 0.9° by 1.25° resolution model configuration. Since the 
resolutions of the NR and the simulated observations are much finer than the CAM- 
Chem resolution, we will use the reduced resolution NR (0.5°, i.e., 42 km 
approximately). The reduced NR simulations are the same as the native NR 
simulations, but the horizontal resolution has been reduced a posteriori (see Da 
Silva et al., 2014). 

Figure 12 displays the reduced resolution (42 km) simulated observations. 
As explained in section 3.5, because the model resolution is 42 km we assume that 
the pixel size has the same size. To generate an appropriate sampling according to 
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the pixel resolution, we divide by a factor of 5 the number of latitude and longitude 
pixels provided in the table. The left panels show the average surface-700hPa 
retrieved CO column for July 2006. The right panels show the cloud free ratio for 
July 2006. For the cloud fraction threshold, we use 20% to keep the same cloud free 
ratio as for the high-resolution observation simulation, as explained in section 3.5. 
Cloud free ratio maps (figure 12) at low resolution are then very similar to the same 
maps at high resolution (figure 6). 

6. Conclusion 

This paper is Part 1 of a two-part study. Here, we demonstrate the feasibility 
of simulating a GEO constellation for air quality monitoring, with a focus on CO. 
Three potential instruments are simulated covering the three most populated and 
polluted areas of the world: Continental US (CONUS), Western Europe and Easter 
Asia. We use very high-resolution output (0.06°, i.e., ~7 km horizontal resolution) 
from the GEOS-5 model as a NR to simulate a MOPITT-like instrument. Instead of 
using a full radiative transfer model to simulate the instrument vertical retrieval 
sensitivity as defined by the AK, we use a novel method described by Worden et al., 
(2013). This method employs multi-linear regression using predictors (from the 
NR) and training set coefficients (from real MOPITT data) to produce scene- 
dependent AKs, thus allowing a very fast computation of the instrument synthetic 
measurement dataset. Thus, we avoid the computational burden of using a full 
radiative transfer model, allowing the generation of one month of GEO constellation 
data in less than 12 hours. This makes simulation of the GEO constellation 
measurement computationally feasible. The main conclusions of this work are as 
follows: 

1. Instead of using the model resolution as the instrument pixel resolution, and 
the defined field of view as a simple latitude/longitude rectangle, we present 
a method to simulate the data using a GEO projection. This gives accurate 
GEO instrument spatial resolutions and fields of view that vary with latitude 
and longitude. 

2. This paper extends application of the Worden et al., (2013) averaging kernel 
(AK) prediction method. Realistic variations of potential GEO instrument 
vertical retrieval sensitivities are simulated. Instrument sensitivities depend 
on predictors and the main drivers are: surface pressure, CO profile and 
temperature profile. Rather than using an average AK for fast computation, 
the observation simulator presented here is able to provide fast computation 
of AK variability (and its associated retrieval error covariance matrix) at the 
same time. 

3. We discuss limitations of the method used for this study. The very high CO 
concentrations occurring in the NR over very polluted areas often overreach 
the training set statistical coverage. In this situation, we use a tuning method 
to reduce the amplitude of CO variations in the NR. 

4. To make the observation simulations as realistic as possible, we account for 
the impact of clouds. Cloud contamination in the observations is strongly 
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dependent on the instrument spatial resolution and the geographical region 
of interest. The Mediterranean weather regimes show the lowest cloud 
occurrences, whereas subtropical weather might provide comparatively 
lower temporal and spatial sampling for air quality GEO measurements. 

5. We present case studies for the three measurement domains and show that 
the observation simulation method employed here provides realistic AK 
variability. The degrees of freedom for signal for the lowermost troposphere 
(DFS 0 - 3 ) ranges from 0.2 to 0.7 with significantly larger values over land and 
for day that reflect the enhanced vertical sensitivity possible with 
multispectral retrievals. We simulate small local DFS 0-3 variations according 
to surface and atmospheric parameters (e.g., surface pressure, CO profile and 
temperature profile). 

6. Simulated retrieval errors that are derived from the AK simulation are 
compared to the true smoothing error. Comparisons show that the retrieved 
errors are realistic, being lower than or in the range of the smoothing error. 

The next step in this study (Part 2) will be to assimilate the synthetic 
measurement data into a global model. To do so, we present here an additional set of 
simulated observations at a reduced spatial resolution (42 km). This allows an OSSE 
for the potential future prediction system of global air quality with the same 
capabilities for each region of interest: the same models (NR and CR), the same data 
assimilation system (AS) and the same instrument design (observation simulator). 
The goals of Part 2 will be to: (1) assess the ability of the GEO constellation to observe 
the impact of pollutant emissions over each region; (2) look at the importance of 
long-range transport between regions; and (3) investigate the value of the 
measurements from each mission in the GEO constellation, taken individually and 
together. 
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874 Tables: 

875 

876 

877 



GEO-AM 

GEO-EU 

GEO-AS 

Subjon 

-97° 

CO 

4^ 

o 

120° 

Number x pixels 

500 

400 

400 

Number y pixels 

230 

250 

200 

Xmax 

3.5° 

2.4° 

3.3° 

Xmin 

-3.5° 

-2.4° 

-3.3° 

Ymax 

7.2° 

8.2° 

6.7° 

ymin 

4.2° 

5.7° 

3.5° 


878 Table 1. GEO-constellation instrument specifications: satellite position, 

879 number of pixels and angles of views. 

880 
881 
882 

883 

884 

885 

886 


Predictors 

Training sets 

9sza 

North Hemisphere Ocean (TIR) 

Emissivity 

CONUS Day (Psrf>900hPa, TIR+NIR) 

Latitude 

CONUS Night (Psrf>900hPa, TIR) 

Surface temperature 

Europe Day (Psrf>900hPa, TIR+NIR) 

dP = P surface" P ref 

Europe Night (Psrf>900hPa, TIR) 

CO column 

Eastern Asia Day (Psrf>900hPa, TIR+NIR) 

Water Vapor Q(z) 

Eastern Asia Night (Psrf>900hPa, TIR) 

COfz) 

N.H. Mountains Day (900hPa>Psrf>800hPa, TIR+NIR) 

Thermal contrast dT(z)=(Tsrf-T(z)) 

N. H. Mountains Day (800hPa>Psrf>700hPa, TIR) 

dT(z) 2 

N.H. Mountains Night (900hPa>Psrf>800hPa, TIR+NIR) 

Cos(0 S za)/logioCO(z) 

N. H. Mountains Night (800hPa>Psrf>700hPa, TIR) 

dT(z) )/logioCO(z) 



88/ *Pref=1000hPa 
888 

889 Table 2. Left: List of predictors, right: List of the different training sets used to 

890 produce the geostationary constellation CO measurements. TIR and NIR state 

891 if the training set simulates multispectral or TIR-only retrievals (see text for 

892 details). 

893 

894 

895 

896 

897 

898 

899 


20 




900 

901 

902 

903 

904 

905 

906 

907 

908 

909 

910 

911 

912 

913 

914 


ACCEPTED MANUSCRIPT 


Figures: 



Figure 1. The chemical OSSE framework. See text for details. 
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Figure 2. Total carbon monoxide 
column in molecules/cm 2 from the 
GEOS- 5 7km resolution Nature Run, 
15 July 2006 15:00 UT. 
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916 

917 Figure 3. Surface CO time average during July 2006 over (a) North America, (b) 

918 Europe, and (c) Asia. 
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922 Figure 4. a) Geostationary constellation measurement domain a) Polar 

923 projection, b) GEO-EU domain in a geostationary projection, red dots are the 

924 full resolution footprints, purple dots are plotted every 100 th pixels, c) is the 

925 same as b) but in an equidistant latitude-longitude cylindrical projection, d) 

926 Geometrical sketch of the geostationary projection. 

927 
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a) MOPITT AVK (DFS= 1.91589) b) MOPITT AVK (DFS= 1.50446) 




928 

929 Figure 5. Examples of original typical MOPITT averaging kernels (AKs). Left 

930 panel: multispectral day/land AK. Right panel: night/land or sea AK. 
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Figure 6. Cloud free ratio (%) for the three measurement domains during July 

2006. 
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Figure 7. Sensitivity matrices of the average cloud free ratio (in %) for pixel 
size versus cloud fraction threshold. 


a) Pixel size 7km, Cloud Fraction Treshold 0.05 


b) Pixel size 7km, Cloud Fraction Treshold 0.25 



939 

940 

941 Figure 8. Examples of cloud detection and ratio of observed area for two 

942 different cloud fraction thresholds and two different pixel sizes. Red are cloud 

943 contaminated pixels and blue are cloud free pixels. Performed over South East 

944 CONUS 5 July 2006 00UT. 
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Figure 9. Snapshots of the Nature Run surface to 700 hPa partial column (a, d, 
g). Corresponding retrieved partial column (b, e, h) and corresponding 
degrees of freedom for signal (DFS) for surface to 700 hPa (c, f, j). Snapshots 
are captured at daytime but different dates following regions: 4 July 2006 
02UT CONUS, 14 July 2006 10UT Europe, 22 July 2006 18UT Eastern Asia. 
Deep colors are the cloud-free pixels. Faded colors represent cloud- 
contaminated pixels that are not used in further processing. 


28 





ACCEPTED MANUSCRIPT 




0.3 0.5 1.0 2.0 3.0 

CO (10~18 molec/cm2) 


‘ % . 


-60 40 -20 0 20 40 

dP (hPa) 


-10 -5 0 5 10 15 20 30 

dT (K) 


0.7 

0.6 

to 

i 0.5 

8 0.4 
ro 

g 0.3 
Q 

0.2 

0.1 

0.7 

0.6 

<u 

o 0.5 

^ 0.4 

ro 

o 

to 0.3 

LL_ 

Q 

0.2 
0.1 

0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 


957 

958 Figure 10. Scatter plots showing variation of degrees of freedom for signal of 

959 surface to 700 hPa versus predictors with highest impacts to the multi-linear 

960 regression fit (see table 2 and text for details). Red are day-time values (3pm 

961 local time) and blue are night-time values (3am local time) 5 July 2006. 

962 Vertical solid lines indicate the mean value of the distribution used to build 

963 the training sets and dashed lines indicate associated ±a (standard deviation). 

964 Dotted lines indicate associated +2a for CO training set. 
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a) CONUS 


b) CONUS 






Figure 11. Left panels: scatter plots of Nature run surface-700 hPa partial 
columns (Xt) versus corresponding retrieved partial columns (Xr). Right 
panels: Smoothing error (Xs) versus corresponding retrieved error (Xe). Dates 
are the same as described in figure 10. 
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Figure 12. Low-resolution observation simulations used for the assimilation 
runs. Left panels: July 2006 average retrieved CO surface-700 hPa partial 
column. Right panels: Cloud free ratio for July 2006. 
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Highlights 

• A constellation of geostationary platforms for mapping pollutant sources 
and variability is described 

• Observation simulation without radiative transfer model is 
computationally cheap 

• Impacts of clouds are diagnosed and is dependent of the weather regime 

• A detailed analysis of the simulated observation sensitivity is performed 

• Limitations of the method are discussed 




